Hausdorff dimension of critical fluctuations in abelian gauge theories 
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The geometric properties of the critical fluctuations in abelian gauge theories such as the Ginzburg- 
Landau model are analyzed in zero background fleld. Using a dual description, we obtain scaling 
relations between exponents of geometric and thermodynamic nature. In particular we connect the 
anomalous scaling dimension rj of the dual matter field to the Hausdorfi' dimension Dh of the critical 
fluctuations, which are fractal objects. The connection between the values of 77 and Dh, and the 
possibility of having a thermodynamic transition in finite background field, is discussed. 
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Anderson has proposed the breakdown of a general- 
ized rigidity associated with prohferation of defect struc- 
tures in an order parameterj-as a general means of char- 
acterizing phase transitionall. In the context of three- 
dimensional superfluids and extreme type-II supercon- 
ductors, such, ideas have recently been put on a quanti- 
tative leveltla. It has been explicitly demonstrated that 
in three spatial dimensions abelian gauge theories such 
as the Ginzburg-Landau theory describing type-II super- 
conductors, suffer a continuous phase transition driven 
by a proliferation of topological defects in the order pa- 
rameter, which are closed loops of quantized vorticityu. 
These loops are induced by transverse phase fluctuations 
in a complex scalar order parameter. Such fluctuations 
are proniincDt in, for instance, doped Mott-Hubbard 
insulatoraj'clS. 

In this paper, we investigate the non-trivial geometric 
properties of these critical fluctuations, and give a geo- 
metric interpretation of the anomalous scaling dimension 
of the condensate order parameter both for a charged and 
neutral condensate. In addition, we discuss the connec- 
tion between the geometric properties of the zero-field 
critical fluctuations and the possibility of having a ther- 
modynamic finite- field phase transition involving unbind- 
ing of loops of quantized vorticity. 

We emphasize that the main results to be presented 
are quite general, and apply to the static critical sector 
of any theory of a complex scalar matter field coupled, to 
a fluctuating gauge-field in three spatial dimensionsav'U, 
provided the symmetry group of the theory is abelian. 

The Hamiltonian for the system is given by 

1 



H{q, u^) = ml + ^ 101^ + iD^^r + -F' 



(1) 



where = F^.F^", F^, = d^h,~ 

and (f) = 101 exp(i0) is a complex matter field coupled to a 
massless gauge field h with coupHng constant q. The |</)|*- 
term mediates a short-range repulsion, while the gauge- 
field h mediates long range interactions, m^ is the mass- 
parameter for the 0-field, and is a self coupling. 

Consider Eq. |l| representing a SD condensate with 
charge q ^ sustaining stable topological objects in the 



form of closed vortex loops. Then the theory with q = 
is a field-theoretical description of the ensemble of these 
stable topological objects, constituting a dual description 
of the original theorya. The theory with g = is also a di- 
rect field-theoretical description of a neutral condensate. 
Thus, in 3D, the gauge-theory H{q ^ 0,u^) describing 
a charged condensate has field-theoretical description of 
its critical fluctuations or topological defects in terms of 
a theory isomorphic to H{q — 0, u^) describing a similar 
but neutral condensate, and vice versaQ. In this sense, a 
charged condensate has neutral vortices with only short- 
ranged steric interactions, while a neutral condensate has 
charged vortices with long-ranged interactions. In the 
former case, the long-ranged interactions between vortex 
segments are rendered short-ranged by fluctuations of the 
gauge-field in the original theory, i.e. the dual gauge-field 
is massi'sse with mass given by the charge of the original 
problems. 

The anomalous dimension rj for the field is defined 
via the relation 
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where Q{z) is some scaling function, ^ is a correlation 
length, and d is the spatial dimension of the system. This 
correlation function has a geometric interpretation, yield- 
ing the probability amplitude of finding any particle-path 
connecting x and y. In the present work, the particle tra- 
jectories correspond to vortex lines. 

For a random walk of length N iu,d = 3, the probabil- 
ity of going from x to y is given byi3 



P(x,y;iV) 
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The correlation function G(x, y) of the corresponding 
gaussian field theory is found by summing up P(x, y; N) 
for all N 



G(x,y) =^P(x,y;7V) 
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Comparing this with Eq. we find 77 = 0, as expected. 
The random walker traces out a fractal path with Haus- 
dorff dimension Dh- Moreover, in general the distance 
between two points x and y N walks apart is given by 

(|x-y|2)cxiV2^, (5) 

where A is the wandering exponent which for the gaus- 
sian 3-D case is A — 1/2. Inverting Eq. ||, we find that 
the total length of the random walker scales with linear 
size as L^/'^, hence the Hausdorff dimension of the ran- 
dom walker is given by Dh = 1/A. If we set x = y in 
Eq. ^ we find that the unnormalized distribution D{N) 
of loops of perimeter iV, at the critical point, is given by 

X 

with a — 5/2 for purely random walkers. The extra fac- 
tor TV^i in Eq. | comes from the arbitrariness in defining 
the starting position along the loop. Hence, for the case 
of strict random walkers in 3D, described by a gaussian 
field theory H = m^|(/)p -I- |V(/)p, the corresponding set 
of values for the two geometric and one thermodynamic 
exponents is given by (A, a, 77) = (1/2,5/2,0). 

Beyond the gaussian case, exact exponents can not be 
obtained analytically, however we will derive scaling re- 
lations for them. When Eq. || is invoked, a generalized 
probability function P(x, y; N) may be written on the 
form 

where F(x) is a scaling function, and normalizability of 
P implies p = dA. From Eq. we find that P(x,y;iV) 
should scale with N as N^~°', which yields the scaling 
relation p — a ~ 1. Conversely, summing over all N in 
Eq. 1^ to find the correlation function G{x, y), we obtain 

G(x,y) =^F(x,y;7V) cx (8) 

N |x-y|^ 

giving the scaling relation 77 = + 2 — rf. Combining 
the above, we find 

rj + DH = 2, Dh = ^—. (9) 

a — 1 

A computation of the geometric exponent a yields the 
thermodynamic exponent 77 and the Hausdorff dimension 
Dh- Note that both 77 and Dh are sensitive functions 
of a, dri/da = —dDn/da ~ d/{a — 1)^, such that a 
precise determination of 77 and Dh requires great pre- 
cision in the determination of a. The above reinforces 
the statement that a geometric transition of the vortex 
tangle at criticality of the gauge theory Eq. |l| can be 
assigned a genuine thermodynamic order parameter via 
a dual formulation of the original theory in three spatial 



dimensionsa^u. The random walker is represented by a 
gaussian theory, Eq. |l| with {u^ — 0, q — 0), for which 
7; = 0. This corresponds to Dh — 2, such that the ran- 
dom walker in three dimensions traces out a path that 
precisely fills a cross-sectional area of the system. Note 
that Dh < 2 ^ 77 > 0, while Dh > 2 ^ 77 < 0. 

The Hamiltonian Eq. |l| with {u^_^ 0, (7 ^ 0) has a dual 
field theory corresponding to Eq. |l| with q — Q describing 
the neutral vortex tangle of a charged superconductom. 
The l0|*-term in Eq. [^represents a steric repulsion^ i.e. 
the vortex loops can not overlap, leading to a random 
walk problem with self-avoiding links (but not necessar- 
ily self-avoiding sites), in the sense that parallel vor- 
teX|-segments repel, perpendicular vortex segments can 
cutcl, while antiparallel vortex-segments can annihilate. 
Hence, this is not a standard self-avoiding path problem. 
However, we expect A > 1/2 or equivalently Dh < 2, 
since steric repulsion should result in a vortex- loop tangle 
packing space less densely than for the non-interacting 
case, so that 77 > 0. The repulsive interaction between 
parallel vortex segments also leads to a more efficient sup- 
pression of long loops than for the non-interacting case, 
so that a > 5/2. 

Consider next Eq. |^ with {u^ ^ 0,(? = 0) for d = 3, 
which has a dual field theory corresponding to Eq. |l| 
with q ^ describing the charged vortex tangle of a 
neutral condcnsateO. A long-ranged (anti) Biot-Savart 
interaction is mediated by the gauge-field. This is a rel- 
evant perturbation, in renormalization group sense, to 
a steric contact repulsionEI. The geometric properties of 
the charged vortex tangle are a result of a balance be- 
tween attractive forces mediated by the gauge field, and 
the steric repulsion. As the numerical simulations show, 
we find A < 1/2, corresponding to Dh > 2 which means 
that the vortex tangle is more compact than the ensemble 
of pure random walkers, due to the fact that an attractive 
long-ranged Biot-Savart interaction between oppositely 
oriented vortex segments overcompensates the steric re- 
pulsion so as to contract the vortex-loop tangle not only 
compared to the pure |0|'*-case, but even compared to the 
noninteracting case. The tangle thus packs space so that 
it more than fills a cross-sectional area of the system. 

The fiuctuation-dissipation theorem provides a bound 
on 7/ via the susceptibility X4> = J d'^xG{x) ^ , which 
is bounded by the volume L"^ of the system, i^^^' — 
j^d _ j^2-d-r) ^ j^d^ so 77 > 2 — d. Eq. ^ gives a geometric 
interpretation of this bound. Specializing to d = 3, 77 = 
— 1 corresponds to topological excitations with Dh = 3, 
an upper limit. 

For d — the continuous phase transition in a su- 
perfluid or extreme type-II superconductor has recently 
been derrumstrated to be driven by a proliferation of vor- 
tex loopsaB. From the above, 7; = — 1 means that a single 
vortex loop at Tc packs space completely, i.e. its perime- 
ter N scales as iV cx L"^, implying that the vortex-tangle 
collapses on itself, rendering the transition discontinuous. 
This may be been seen from the standard scaling rela- 
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tion (3 = V {d—2 + r])/2 for critical exponents. Formally, 
this implies that the limit r] — > (2 — d)"*" corresponds 
to the limit /3 ^ 0+, characteristic of a discontinuous 
transition. More informally, a collapse of a vortex tan- 
gle may be viewed as mediated by an effective attractive 
vortex interaction, a situation akin to what is known in 
type-I superconductors. Deep in the type-I regime, it is 
known that svioerconductors suffer a weakly discontinu- 
ous transitionll3. 

Monte Carlo simulations have been performed on the 
lattice version of Eq. |l| in the phase-only approximation, 
to determine precise values of a, both for g = and q ^ 
0. We have also performed simulations on pure random 
walkers described by the theory H{q = 0, = 0). They 
reveal that a determination of a is less fraught with finite- 
size effects than a determination of Dh- Thus, we focus 
on determining a. The model we consider is 

H^-J cos(0. - - qc h,j) + i(V X h)2, (10) 

<ij> 

where the site-variable 6i is the phase of the complex 
matter field = |0|exp(z6') of Eq. |l|, when the system 
is discretized, J is essentially a bare phase-stiffness, and 
the link- variable hij = d\ ■ h.. The charge qc is the 
(original) charge entering in the simulations. Up to this 
point we have considered a general charge q irrespective 
of whether it couples to the original condensate or the re- 
sulting vortex tangle. The numerical simulations are per- 
formed on the phase of the condensate, hence the concept 
of original and dual are fixed in terms of the numerical 
simulations. Consequently we introduce the charges gc 
for the condensate and qy for the vortices. 

^From the phase distrilDutions of the matter field we 
can extract vortex loopsa. These loops have charge qw 
and are described by the field theory H{qy, Utpp. Hence, 
we can study the critical properties of the charged field 
theory H{qY, u^f,) by considering the geometric properties 
of the thermally excited vortex-loop tangle at the criti- 
cal temperature in the 3DXY model. Conversely, the 
geometric properties of the vortex tangle with gc 7^ 
yield the critical properties of the neutral field theory 
if(qv = 0, u^) of Eq.J. n 

The simulations with qc — ^ are described elsewhereB, 
while for qc ^ ^ the simulations proceed as follows. 
For every site on the lattice a phase change 6i 9[ 
is attempted, and accepted or rejected according to the 
Metropolis algorithm. Then a change in hij — > hij + Sh 
is attempted, and accepted or rejected according to the 
Metropolis algorithm. When updating hij we update all 
the link variables on a randomly oriented elementary pla- 
quette containing hij as one of its four edges. Updating 
of h in this fashion guarantees that the gauge-fixing con- 
dition V ■ h = is enforced at all times. For qc = 0, the 
simulations were performed for a system of size Lx Lx L 
with L = 180, while those for qc ^ were performed with 
L = 64. 



During the simulations we have sampled the distribu- 
tion function D{N), Eq. ??, obtaining a. The results 
are shown in Fig. |l] and listed in Table I. The value of a 
obtained for qc 7^ (dual neutral), which is the hardest 
system to simulate, gives a value for 77 in good agreemeat 
with high-precision results for ry of the pure 0'*-theorylliJ. 
This serves as a useful benchmark on our method of ex- 
tracting T]. For qc = we have simulated much larger 
systems than for qc 7^ 0. The deviation from the gaus- 
sian value a = 5/2 is substantial, and of opposite sign 
compared to qc 0. Given the size of the system we 
consider for qc — 0, it is unlikely that this is a finite-size 
artifact. An a < 5/2 guarantees 77 < for the qc — 
(dual charged) case, contrary to the value of 77 > for 
qc 7^ OIIj'cI. In particular, the inset of Fig. 1 lends strong 
support to the proposition that ri{qc 7^ 0) > 0, while 
7]{qc = 0) < 0. 

The value 77 < obtained for the original neutral, dual 
charged case, is significant: It implies that Dh > 2 for 
this case. Whether Dh > 2 or Dh < 2 is of great import 
to the possibility of having a genuine phase-transition 
driven by a vortex-loop unbinding even in the presence 
of a finite background field such as magnetic induction in 
type-II superconductors. A vortex system accesses config- 
urational entropy more easily if it is compressible than if 
it is incompressible. For the charged case the gauge-field 
fluctuations render the system compressible. In the neu- 
tral case, the system expands screening strings of closed 
vortex loops to a larger extent than for the charged case, 
as substitutes for the gauge-field fluctuations. This is 
why Dniqc = 0) > Dniqc 7^ 0). There is an infinitely 
larger amount of screening vortex-strings in the neutral 
case (dual charged) than for the charged case (dual neu- 
tral), which is the true significance of the fact that rj 
is smaller for qc = Q than for qc ^ Q. The possibility 
of the zero-field vortex-loop blowout transition surviving 
the presence of a finite field is much greater in a neutral 
superfluid or an extreme type-II superconductor, than in 
a charged condensate with a priori good screening. 

Given the significance of Dh > 2, we elaborate on the 
fact that for the neutral (dual charged) case, we find -q < 
0. The Lehmann-representation of the Fourier transform 
G{p) of Eq. H, is sometimes used to argue that ij obeys 
the strict inequality 77 > 0. The Lehmann-representation 
of G{p) is given by 

G{k)^ Td^,' (11) 
Jo -\- 

where 1 = dfi^ pi^ji^), and p{^P') = Z5{p? — m^) -I- 
(t(/i^). The propagator for the gaussian case would be 
G{k) — l/(fc^-|-?Ti|), where ?7i^ refers to the bare masspa- 
rameter in Eq. ^ Thus, 77 > follows if < Z < 1, which 
holds for a uniformly positive p(/^^). However, in theo- 
ries with a local gauge symmetry, the two-point correla- 
tion function is a gauge-dependent quantity. Thus, cr(//^) 
may in principle be made negative for certain values of fj, 
by a gauge-transformation. This invalidates the reason- 
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ing leading to the strict inequalities Z < 1 and rj >|-0-| 
A negative ij, as found here and in other simulationaj'u 
all representing basically exact results- agrees with a re- 
cent non-perturbative RG calculationllj, which also gives 
I' = i^sDXY at the charged critical point. 

At the critical point, the relevant flucjtjuations arc 
transverse phase-fluctuations, or vorticeaau'B. Ignoring 
amplitude fluctuations yields an effective Hamiltonian 
governing the transverse 0- fluctuations, whose Fourier- 
transform !F we denote by Sfe, J^{{'V9)t) = Sk ~ 
—2tt i (k X n/;)/fc^, where is the Fourier-transform 
of the local vorticity. We find, after integrating out the 
transverse gauge-field, that H = S^(fc) Sk ■ S_fc, where 
S2(fc) = fcV(yfc2 + 2q2). For q = 0, we have E^{k) = f , 
while limfc_,oS^(fc) ~ for q 0. The coupling to 
a fluctuating gauge field softens the transverse phase- 
fluctuations, providing the effective phase-stiffness with 
an extra power compared to the q = 0-case. Thus, 
G-i(fc, q = 0) = k^+^{k) and G-\k, q^O)^ k'^+^{k). 
In both cases, the k 0-limit of the self-energy S(fc) is 
given by S](fc) ~ k^^^ . We thus have \imk-,o (k) ~ 
provided ry > for q — and, when invoking 
the absolute lower bound, rj > —1 for q 0. For 
a pure |(/)|*-theory, the Lehmann-reprcsentation coupled 
with positive-deflniteness of p{is^), holds. 
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FIG. 1. The vortex-loop distribution function 

D{N) ~ N^" as a function of loop-perimeter A'^, at the 
critical point for the charged (qc 7^ 0, qv = 0) and neutral 
(gc = 0, gv 7^ 0) cases. Numerical results for the exponents 
{a, Dh, a, rf) are given in Table I. The system size is Lx Lx L, 
with L — 180 for qc ~ 0, and L = 64 for qc 7^ 0. Inset 
shows simulation results for N^^'^ D{N) ~ ]\[^/^-'^ on a dou- 
ble-logarithmic scale. Top: qc ~ 0, L = 180 (dual charged). 
Middle: Noninteracting vortex loops (gaussian case) , L = 64. 
Bottom: qc 7^ 0, L = 64 (dual neutral). The results demon- 
strate that 5/2 — Q > for qc ~ (dual charged), while 
5/2 — Q < for qc 7^ (dual neutral). Hence, by Eq. ^, 
Dh > 2,ri < for qc = 0, while Dh < 2,ri > for qc / 0. 
The latter agrees with other high- precision results for 77, see 
Ref. 11. Note that the gaussian result a = 5/2 is obtained to 
high precision, for L — 64. 



^ P. W. Anderson, Basic Notions in Condensed Matter 
Physics, Benjamin Cummings, (1984). See also P. G. de 
Gennes, Proceedings of the Nobel Symposium, 24, 112 
(1973). 

^ Z. Tesanovic, Phys. Rev. B 59, 6449 (1999). 

^A. K. Nguyen and A. Sudb0, Phys. Rev. B 60, 15307 

(1999); Europhys. Lett., 46, 780 (1999). 

V. J. Emery and S. A. Kivelson, Nature, 374, 434 (1995). 

See also H. Kleinert, Phys. Rev. Lett., 84, 286 (2000). 
^ J. Hove and A. Sudb0, Phys. Rev. Lett., 84, 3426 (2000). 
® H. Kleinert, Gauge fields in condensed matter, Vol. 1, 

World Scientific, Singapore, (1989). 

D. Austin et al, Phys. Rev. D 49, 4089, (1994). 
* A. Sudb0 and E. H. Brandt, Phys. Rev. Lett. 67, 3176 

(1991); C. Carraro and D. S. Fisher, Phys. Rev. B 51, 534 

(1995) . 

^ C. Dasgupta and B. I. Halperin, Phys. Rev. Lett., 47, 1556, 
(1981). 

B. I. Halperin et al, Phys. Rev. Lett., 32, 292 (1974). 
" M. Hasenbusch and T. Torok, J. Phys. A 32, 6361 (1999). 
" I. F. Herbut and Z. Tesanovic, Phys. Rev. Lett. 76, 4588 

(1996) . 



Exponent 


Gaussian 


9C = 0, qv 7^ 


gc / 0, qv = 


Q 


5/2 


2.312 ± 0.003 


2.56 ± 0.03 


Dh 


2 


2.287 ± 0.004 


1.92 ± 0.04 


A 


1/2 


0.437 ± 0.001 


0.52 ± 0.01 


V 





-0.287 ± 0.004 


0.08 ± 0.04 



TABLE I. The loop distribution exponent a, as deter- 
mined from Monte-Carlo simulations. The remaining expo- 
nents have been determined from Eq. ^ Symbols are ex- 
plained in the text. 
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